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ABSTRACT 

When Cold Dark Matter (CDM) power spectrum normalized by COBE results, its 
amplitude at smaller scales can be parametrized by T ~ Qgh. The expected variations 
of Hubble flows in two samples, the sample of 36 clusters in the Mark III catalogue, 
and the sample of 20 Type la supernovae (SNe), are calculated for the power spectrum 
in critical-density CDM models (including tilted CDM models and vacuum-dominated 
CDM models). The comparison between the expectations and the real variations in 
the data offers a bias-free way to constrain the power spectrum. The cluster sample 
yields T < 0.30 - 0.88(n ps - 1) + 1.9(ra ps - l) 2 at 95% C.L., with best fits T = 
0.15 — 0.39(n ps — 1) + 0.37(rt ps — l) 2 , where n ps is the spectral index of the power 
spectrum. The Type la SN sample yields T < 0.25 - 0.80(n ps - 1) + 1.6(n ps - l) 2 
at 95% C.L., strongly favoring lower T's. The results are inconsistent with a critical- 
density matter-dominated universe with n ps ^ 0.8 and Hq J; 50 km/sec/Mpc. 
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1 INTRODUCTION 

Variations in Hubble flows (peculiar Hubble flows, or 
monopolar deviations from a global Hubble flow) are directly 
tied to the underlying density fluctuation. For example, for a 
sphere with a uniform density embedded in a homogeneous 
FRW universe with a critical density, its peculiar Hubble 
flow under linear theory is 



ACDM Models where A is the cosmological constant) is 
parametrized by 



P(k) = 2Tr 2 (5 2 (3000/i _1 Mpc) 3 



(3) 



5H_ 

H 



ISp 

37' 



(1) 



where SH/Ho is the local deviation from a global Hubble 
flow normalized by the Hubble constant Ho, and 5p/p is 
the local density contrast. For realistic Slo = 1 models with 
gaussian power spectra, the r.m.s. peculiar Hubble flow as 
measured within a top-hat sphere is, however (Turner, Cen 
and Ostriker 1992; Shi, Widrow and Dursi 1995) 



where (Bunn and White 1997) 
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and (Bardeen et al. 1986; Sugiyama 1995) 
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Thus, if we know the variations of Hubble flows as a func- 
tion of scales, we have a grasp of the power spectrum that 
generates the underlying density fluctuations. Since only pe- 
culiar velocities are involved, the approach is free of the bias 
factor b. And ideally, it will not be influenced by structures 
beyond the sample volume. 

Conventionally, the power spectrum P(k) in critical- 
density CDM models (including tilted CDM models and 



In the equations, flo is the matter content of the universe, 
fii, (adopted to be 0.024h" 2 , Tytler, Fan and Buries 1996) 
is the baryonic content of the universe, and h is the Hub- 
ble constant in the unit of 100 km/sec/Mpc. Eq. (Q) is due 
to normalization by the COBE results. Spectral index n ps 
is constrained to between w 0.8 and 1.1 (Lineweaver and 
Barbosa 1997). Given the COBE normalization, T and n ps 
determine the shape of the power spectrum. 

In this article we calculate the expected variations of 
Hubble flows for real samples using the above form of the 
CDM power spectrum. We then compare the expectations 
to real variations and obtain constraints on T as a function 
of n DS . 
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2 FORMALISM 

A detailed formalism can be found in Shi (1997). Here we 
only present some essential equations. 

The deviation from a global Hubble flow, 5H (= 
^Locai — Ho), of a sample and its bulk motion U are 



6H = B~ 1 sp S q r q -U irq r\ 
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where 

q q q' <i 

Vector v q is the position of object q in the sample (with us 
at the origin), and S q is its estimated line-of-sight peculiar 
velocity with an uncertainty a q . 

The expectation of the deviation 5H consists of two 
parts: the contribution from the density fluctuations 5H^ V \ 
and the contribution from noises 

In linear theory, 
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where the window function of the sample 
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Spatial indices i,j,l,m run from 1 to 3. Identical indices 
denote summation. 

The r.m.s. noise contribution is 
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Since we do not know the precise value of Ho, we can 
only deal with relative variations of Hubble flows within a 
sample. If we denote the expansion rate defined by a sub- 
sample with n nearest objects (relative to us) as H n , that 
by another subsample with m(< n) nearest objects as H m , 
and that by the entire sample with N objects as Hn, then 
the expected 

H n — H m \ 2 \ 
H N J I 



in which 
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and 
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To better account for the contribution from fluctua- 
tions on very small scales at which the window functions are 
choppy and random (see figure 1) and density fluctuations 
become non-linear, we use the variation of 5H as a measure 
of the amplitude of density fluctuations, and represent the 
small-scale contribution by a one-dimension gaussian ran- 
dom motion with a r.m.s. velocity cr». The expectation is 
then 
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where 5 nm is the Kroeneker 5 function. Thus, given AH„ — 
5H n — 8H n -\ from a real sample D, the probability of a 
critical-density CDM model with a particular set of Y and 
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3 RESULTS 

The formalism is applied to the cluster sample that yields 
the template Tully-Fisher relation in the Mark III catalogue 
(Willick et al. 1997). Only subsamples with no less than 15 
clusters are included in calculations, because smaller sub- 
samples are susceptible to large effects of non-linearities. 
This cut-off does not significantly affect the final results due 
to the large uncertainties of SH in the small subsamples. 
The integrals over the window functions are done by sam- 
pling about 40000 k's in space. The window functions are 
truncated at n/k w 6ft -1 Mpc. To include the effect of non- 
linearities below the truncation scale, a a, = 300 km/sec is 
added quadratically to the uncertainty in distance of each 
cluster to fully represent the uncertainty of its radial pecu- 
liar velocity. The truncation scale is verified not to signif- 
icantly influence our final result, but the magnitude of a* 
does have an appreciable impact. A smaller <j» will increase 
the contrast of the probability distribution P(M\DI) , thus 
tighten the resultant constraints. At the moment, a cr» of 
300 km/sec is quite a safe assumption to make (Bahcall and 
Oh 1996). 

The test is also applied to the sample of 20 Type la SNe 
of Riess, Press and Kirshner (1996). Subsamples run from 
the 6 nearest Type la SNe to the full sample of 20 SNe. 
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Truncation of window functions is made at %/k ~ 4h~ 
Mpc. The host galaxies of the SNe are assumed to have a 
<7* = 480 km/sec. This corresponds to a r.m.s. gaussian ran- 
dom motion of 830 km/sec for host galaxies, which roughly 
matches the velocity dispersion of galaxies in clusters. It is 
also the maximal <r* allowed by the sample, because any 
addition of distance uncertainties will result in a dispersion 
of radial peculiar velocities that is larger than that of data. 
Therefore, the constraint we obtain will be rather conserva- 
tive. 

Figure 2 shows the probability distribution (with an 
arbitra ry normalization) of V when Oh — 1 and 



for the I two samples After the same distribution function 



is calculated for other Qo (in which case A = 1 — Qo) and 
n ps , we find that at 95% C.L., the cluster sample yields 

- I) 2 , with best fits T = 
the Type la SN sample 



T < 0.30 - 0.88(n ps - I) + 1.9(n ps 



-If 

+ 1.6(n ps - l) 2 at 95% C.L. 



0.15 - 0.39(n ps - 1) + 0.37(n pa 
yields F < 0.25 - 0.80(n ps - 1) 
strongly favoring lower F's. These results are accurate for 
critical-density CDM models in the range 0.3 < fio < 1 to ~ 
10%, because the factor Oq' 6 in peculiar velocities is mostly 
offset by the fio dependence in the COBE normalization 
eq. (^). The results indicate that a critical-density matter- 
dominated universe with Ho <; 50 km/sec/Mpc and n ps J; 
0.8 cannot fit the variations of Hubble flows in the data. 
(The constraints will not change significantly if one includes 
the tensor contribution to the COBE normalization of P(k), 
but will tend to relax when one decreases P(k) on ~ 100/i -1 
Mpc scales.) With the number of Type la SNe with distance 
measurements doubling very soon, we can only expect that 
the limit on F becomes more secure and interesting. 
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